Inter-filament Attractions Narrow the Length Distribution of Actin Filaments 
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We show that the exponential length distribution that is typical of actin filaments under phys- 
iological conditions dramatically narrows in the presence of (i) crosslinker proteins (ii) polyvalent 
counterions or (iii) depletion mediated attractions. A simple theoretical model shows that in equi- 
librium, short-range attractions enhance the tendency of filaments to align parallel to each other, 
eventually leading to an increase in the average filament length and a decrease in the relative width 
of the distribution of filament lengths. 

PACS numbers: xx.xx 
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Introduction 

The protein actin in its filamentous form (F-actin) is a 
major structural component in the cytoskeleton network, 
and plays an active role in maintaining the leading edge of 
moving cells [l|. In all these structures, the filaments are 
dynamic objects that continuously polymerize and de- 
polymerize as part of their normal function. The physical 
properties of the overall structure (mechanical strength, 
viscoelastic response, etc.) depend on the properties of 
the filaments that form the structure and, in particular, 
on their length distribution. 

In vitro experiments on F-actin under physiological 
conditions (with and without capping proteins) reveal a 
wide, exponential distribution of lengths Q, y, |^ . This 
can be attributed theoretically to the fact that the under- 
lying (dominant) stochastic dynamics, namely monomer 
exchange, is an "homogeneous zero range process" (ZRP) 
|5J which generates an exponential steady state distri- 
bution. At the same time, equilibrium theories of self- 
assembly of polymers or micelles predict an exponential 
distribution determined by the competition of the end- 
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FIG. 1: A schematic view of three different attraction agents 
(" linkers" ) and their effect on actin polymerization: (a) crosslinker 
proteins (b) polyvalent counterions, and (c) inert polymers giving 
rise to depletion mediated attractions. 



cap energy and the translational entropy [l^. Indeed, 
nucleating and severing proteins such as gelsolin [g are 
known to modify the average filament length, (l), but 
not the coefficient of variance of the length distribution 
Tct = ai/{l), where ai is the standard deviation of the 
distribution. 

In contrast, we have recently observed that when 
crosslinker proteins are added to an F-actin solution they 
enhance the preference for overlapping filament sections. 
As a result, the shape of the distribution changes and r^r 
is reduced approximately by half 7^ . Similar behavior is 
also observed in the presence of myosin aggregate Q . 

In this letter, we argue that the narrowing of the length 
distribution is a general effect arising from short-range 
attraction between the filaments; chemically crosslinking 
is not necessary to observe the reduction in r^. Specifi- 
cally, we show experimentally that both depletion me- 
diated attractions in a solution containing inert poly- 
mers (PEG), and electrostatic interactions induced by 
multivalent counterions (spermine) result in attractions 
that significantly narrow the distribution. A simple the- 
oretical approach demonstrates that this effect should be 
observed for any mechanism that leads to short-ranged 
attraction. The attractions increase the tendency of fil- 
aments to lie parallel to each other and to grow concur- 
rently. As a result, the shape of the length distribution 
changes dramatically: from a monotonously decreasing 
exponential distribution, to a non-monotonous Gaussian- 
like distribution with a well-defined peak at large fila- 
ment lengths. Generally speaking, exponential distribu- 
tions where the width scales as the average, are typical of 
systems where entropy dominates over interaction ener- 
gies; Gaussian-like distributions, where the width is much 
smaller than the average, are typical of systems where in- 
teraction energies dominate. 

Conflicting experimental measurements for the values 
of the elastic constants of the actin gel have been sug- 
gested to be a result of variations in the length distribu- 
tion of F-actin Q- Furthermore, cells provide a biochem- 
ical environment crowded with macromolecules. Crowd- 
ing leads to depletion mediated attractions that can be 
reproduced in vitro by adding inert polymers such as 
polyethylene glycol (PEG) 10]. We therefore suggest 
that the physical effects described in this study have a 
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FIG. 2: Normalized fluorescence data from an actin depoly- 
merization measurement. The fuU circles denote F-actin 
which was incubated in the presence of lOOTTiAf spermine and 
depolymerized by diluting the solution 10-fold with 0.4X AP- 
buffer. Empty circles denote a control measurement where the 
F-actin solution (without spermine) was diluted with a buffer 
containing lOmM spermine. The dashed curves show the non- 
parametric fits to the data and the straight lines are drawn 
where the intensity of fluorescence reaches half of its initial 
value. Inset: an example of a cumulative distribution func- 
tions of filament lengths which was extracted from the data. 
The depolymerization rate, rj, is normalized as explained in 



role in regulating the properties of actin gels in cells and 
in tissues. 

Experimental materials and methods 

Spermine (4 HCl) and PEG (MW= 6000) were purchased 
from Sigma- Aldrich Ltd. (Rehovot, Israel). Pyrene la- 
belled actin (10% labelled monomers) was purchased in 
lyophilized form from Cytoskeleton Inc. (Denver CO, 
U.S. A) and polymerized at a stock concentration of 
0.15mg/ml in OAX polymerization buffer with O.bmM 
ATP ("AP-buffer") as described in 0. G-actin aliquots 
were quick frozen in liquid N2, and pre-spun at I5OOOO5, 
4°C, for 2 hours upon thawing. After polymerization, 
equilibrated spermine (lOOmM) or PEG {A%w/w) were 
added to the F-actin solution from lOAT stocks (diluted 
in AP-buffer) and an equal volume of buffer was added to 
control F-actin aliquots. The filaments were then incu- 
bated for at least two hours at room temperature before 
performing a depolymerization assay. 

Depolymerization was induced by diluting the F-actin 
solution 10-fold using a cut pipette tip. F-actin aliquots 
incubated with PEG (spermine) were compared to con- 
trol filaments depolymerized in AP-buffer containing 
0.4% PEG (lOmM spermine). Depolymerization was 
monitored with a Spex FluoroLog-3 spectrofluorometer, 
Jobin Yvon Ltd. (London) at 20°C, while stirring very 
gently to avoid filament breaking. An CD 1 ND filter 
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1.5 ± 0.4 [7] 


20% a-actinin 
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1.1 ±0.2 
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1.2 ±0.3 


5.5 ±1.5 
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0.5 ±0.2 


1.2 ±0.5 



TABLE I: Results from depolymerization measurements un- 
der four different sets of conditions. The "incubation" ("dilu- 
tion") column specifies what was added to the 0.4X AP-buffer 
during the incubation (dilution) of the actin filaments. The 
values for the coefficients of variance and the depolymeriza- 
tion half-times given in each case are averaged over 4 — 6 
experiments for each set of conditions and the errors are sta- 
tistical. 



was used to minimize photobleaching. Photobleaching 
was measured separately and corrected for when found 
to be significant. The data was non-parametrically fitted 
to a smooth, monotonically decreasing and convex curve, 
r^ was then derived from the smooth curve as described 
in ^. 

Experimental results 

Fig. 121 presents data from an actin depolymerization as- 
say in the absence and in the presence of spermine ions. 
We analyzed 5 depolymerization experiments after incu- 
bation with spermine ions, 4 spermine control depoly- 
merizations, 6 depolymerization experiments after incu- 
bation with PEG and 6 PEG control depolymerizations 
(as described in Sec. II). The values extracted for r^ 
and Ti/2 are summarized in Table 1. In both cases. 
To- was about 2-fold smaller and T1/2 was about 3-fold 
smaller as compared with the control depolymerization 
measurements. These results are consistent with results 
obtained for crosslinked filaments y\. When the length 
distribution is narrower the depolymerization process is 
more efficient since the dominant portion of the filament 
length distribution depolymerizes concurrently, resulting 
in smaller values of tj^/2- As explained below, narrow 
length distributions are a signature of attractive inter- 
filament interactions, whether they result from depletion 
forces jll| or from multivalent counterions jl2|: physical 
crosslinking is not a necessity. 

The average amplitude of the fluctuations around the 
smooth fit to the depolymerization curve was similar for 
dilution buffer containing lOmM spermine and not con- 
taining spermine. However, the presence of 0.4% PEG 
in the dilution buffer introduced a noise amplitude ap- 
proximately twice as large. The effect of noise is more 
pronounced at the late stages of the measurement re- 
sulting in an apparent longer tail of the distribution of 
lengths, which leads to an artificially larger value of r^- 
(see [Zl). This is reflected in the fact that the value of 
Tct in the well established case of the exponential length 
distribution of pure F-actin was larger than unity when 



extracted from dilution in the presence of PEG. It is im- 
portant to note that despite this effect, an unambiguous 
decrease of both r„ and T1/2 is preserved. This indicates 
that the interaction does indeed reduce the variance of 
the lengths. 

Theory 

Actin polymerization is a dynamic process and most the- 
ories that address filament length distributions 0, llil [lj| 
are based on kinetic equations for the various processes 
involved (e.g., monomer association, dissociation etc.). 
On the other hand, Flory has already noted _13] that 
quite often, kinetic arguments leading to a steady state 
distribution can be replaced by equilibrium arguments. 
In particular, the presence of actin-sevcring proteins (gcl- 
solin, cofilin, ...) ensures that the filament distribution 
follows the equilibrium distribution and does not depend 
on dynamic instabilities. Indeed, the exponential dis- 
tribution that is typical of F-actin is also found in self- 
assembling systems where surfactant molecules form lin- 
ear micelles and equilibrium is maintained through ex- 
change of surfactant molecules with the solution [l5l| . 
The equilibrium approach has the advantage that it pro- 
vides a natural framework to systematically study the 
effect of inter-filament interactions. 

The theoretical starting point is the following free en- 
ergy (per unit volume) of a solution consisting of concen- 
trations, pi, of filaments of lengths / = 1, 2, 3, ...: 

/ = ^ p; [ln(p/Wo) -l + al + b] 
I 

+ ^w{l,l')pipi, - fi^lpi (1) 

Kl' I 

The logarithmic term represents the translational en- 
tropy of the filaments where vq is a monomer volume. 
Here and in the following, all energies are in units of 
the thermal energy ksT. The energy of a single fil- 
ament consists of a term linear in I which is due to 
the interaction between neighboring monomers and an l- 
indcpendent term that is basically the end-cap energy, b, 
of the filaments. The next term is the two-body interac- 
tion term where w(/, /') = JdrdJ7[l— exp(— u/ ;/(r, H))] is 
the second virial coefhcient of a pair of rods of lengths /, /' 
averaged over their mutual separation r and angle 17 ,16j . 
Finally, the lagrange multiplier p, is added to free energy 
to fix the total monomer concentration ^^ I pi = pm- 

The equilibrium length distribution is obtained by 
minimizing the free energy with respect to pi: 



Pi 



po exp ~al — y^ w{l, l')p{l') 



(2) 



where po — exp(— 5)/vo and a = a — p,. In the absence 
of inter- filament interactions, the distribution p\ ' is ex- 
ponential with an average filament length (Z)o = 1/a = 
{pm/ PoY^"^ and To- = 1. Note, that because of monomer 



conservation the value of a does not affect the length dis- 
tribution and only shifts the monomer chemical potential 
by a constant. 

The leading contributions to the virial coefhcicnt due 
to hard-core repulsions are u'hc(^j ^') = {TT/2)ll'd + n{l + 
l')(P where d <^ 1,1' is the filament diameter |l3|. As 
noted by Ben-Shaul and Gelbart [l3|, the first term 
merely shifts the monomer chemical potential by a con- 
stant and has no effect on the distribution. The next term 
is equivalent to a reduction in po ^^nd leads to an increase 
in the average filament length 17] while the distribution 
remains exponential with r„ = 1. 

This is no longer the case when linker-mediated inter- 
filament attractions are introduced |l8| . The new contri- 
bution reads 



w^tS.l') = -ll'5{l~e-^°) 



(3) 



27rd5 <^ 2/1 



- 1 



{l2-ll)[l-C-'] 



where uq < is the linker mediated short-range attrac- 
tion per monomer and 5 is the range of the attraction. 
The first term is the contribution from configurations 
where the filaments cross at large angles while the sec- 
ond term is the contribution from parallel configurations 
with z = {li/d)uo, h ^ nim(l,l'), and I2 = max(/,/'). 
The latter is divided into a term independent of I2 due 
to configurations where the two rods overlap only par- 
tially and configurations where the short filament is fully 
adjacent to the long filament. 

Eqs. I2I3I allow one to determine the length distribu- 
tion {pi}. In order to simplify the calculation we look for 
a perturbative solution around the noninteracting distri- 
bution of the form pi = p\ {\ -\- 5i) where 5i ^ 1. To 



leading order, the solution is then 



5i = exp I —5a I — V] 'w{l, l')pi, 



(0) 



(4) 



where 5a is determined from the conservation constraint 
Hilpf^^i^Q. ^ 

Typical solutions for the normalized length distribu- 
tion pi ~ pil'Yliv Pv s-i^d the accumulated distribution 
Fi = 'Yl,ii=\Vi' for different values of the attraction 
strength, uq, are shown in Fig. |21 The solutions for 
Mo ^ —0.001 generally agree with the measured distri- 
butions (see inset in Fig. 2) except for the short filament 
end, which is beyond the resolution of the experiment. 
As the attraction becomes stronger, the distribution de- 
velops a peak at longer and longer filament lengths and 
with increasing magnitude. Since uq is defined as the in- 
teraction energy per monomer, it is the magnitude of Iuq 
that determines whether the inter-filament interactions 
are strong enough to modify the distribution. 

The effect on the coefficient of variance r^ is non- 
monotonic as shown in Fig. QJ). At first, the ratio in- 
creases above one but soon it reaches a maximum and 
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FIG. 3: Length distributions pi (a) and accumulated distributions 
Fi (b) for different values of the attraction strength |no| (in units 
of kgT). The values used in the calculation are p,„ = 0.1, po = 
exp(— 10), 5 = 1, d = 1/2tt and the cutoff value for the filament 
length is imax = 1000. 
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FIG. 4: Average filament length (Z) and root mean square A/ (a) 
and their ratio r^ (b) as functions of the attraction strength |uo 
within the perturbative approximation. The results of the mean 
field approximation are shown in (c) and (d) where the horizontal 
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Same physical values as in Fig. |3] 



decreases below one. The initial increase in r^ is due to 
the biniodal length distribution at intermediate values of 
uo, and the consequence decrease occurs when the longer 
filaments dominate the length distribution. 

Because of the perturbative nature of the approxima- 
tion, it can not correctly describe the behavior in the 
presence of strong attractions. In this limit, of strong 



interactions, the behavior can be studied within a mean 
field-like approximation where the main assumption is 
that the sum in the exponent of Eq. |21 is dominated by 
the "typical" filament length l* — (l). It is helpful to 
expand the attraction term w(Z, I') in powers of uq and 
neglect a term proportional to If that is only significant 
in the immediate vicinity of /*. The length distribution 
is then of the form 



Pi 



p^^-al + (9l'p')l^ I < I* 



(5) 



where g ^ ttmq. It is now apparent from Eq.|Slthat while 
the length distribution of long filaments {I > I*) remains 
exponential, that of short filaments (/ < /*) shows two 
local maxima, first at I — and then at I — I* . Fur- 
thermore, while for small values of |uo| the distribution 
reduces to the original exponential one, at large values of 
|uo| the distribution shows narrow peak at / = Z*. With 
three unknowns a,l* ,p* , the distribution function, p/, can 
be calculated by solving the following three equations: (i) 
the monomer conservation condition Jdllpi = pm, (ii) 
the consistency condition p* = poCxp(— a + gl* p*), and 
(iii) the condition that (1) —I*. As shown in Fig. ^ the 
results of this calculation (valid only at large values of 
|uo|) qualitatively confirm the results of the perturbative 
approach. 

Conclusions 

In summary, we have shown that in the presence of short- 
ranged attractions the coefficient of variance of the length 
distribution of F-actin is decreased due to an increasing 
tendency of the filaments to align parallel to each other 
even without chemical crosslinking. The phenomena is 
quite general and does not depend on the exact details 
of the attraction mechanism: the attractions can be in- 
duced by multivalent counterions or can be electrostatic 
by nature; they can occur in a crowded environment like 
the cell as a result of depletion mediated attractions; or 
they can be induced by linker proteins, which attach actin 
filaments to each other in the cytoskeleton and in other 
filament aggregates. 
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